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INTRODUCTION 

In  order  to  obtain  a  comprehensive  model  to  describe  the  interaction  of 
ultrasonic  waves  with  fiber  reinforced  laminated  composite  materials  ,  a  series 
of  experimental  and  theoretical  studies  has  been  undertaken  characterizing  the 
behavior  of  elastic  waves  in  anisotropic  layered  media.  The  principal  motiva¬ 
tion  of  these  studies  has  been  to  develop  improved  nondestructive  evaluation 
techniques  in  which  measurements  are  made  of  properties,  such  as  complex  elastic 
constants,  which  are  critical  parameters  used  in  the  modeling  of  failure  and 
fracture  in  composites . 

This  report  will  be  devoted  entirely  to  the  mathematical  treatment  of  ob¬ 
lique  incidence  wave  propagation  in  layered  anisotropic  media.  Results  will  be 
developed  for  use  in  subsequent  reports  which  will  detail  various  experimental 
studies  and  will  describe  the  use  of  a  comprehensive  computer  code  for  gener¬ 
ating  accurate  numerical  results. 

No  attempt  has  been  made  to  provide  a  comprehensive  treatment  of  wave  prop¬ 
agation  in  all  anisotropic  media,  rather  consideration  has  been  restricted  to 
symmetries  which  are  found  in  certain  composite  materials  of  interest.  Thus, 
the  results  apply  to  materials  of  orthorhomic  symmetry  and  higher.  Within  this 
class  of  problems  a  further  restriction  has  been  made  limiting  boundary  condi¬ 
tions  to  treat  only  those  boundary  planes  coincident  with  principle  planes  of 
symmetry  for  the  material.  In  so  doing  we  eliminate  a  few  cases  of  practical 
interest  but  obtain  the  advantage  of  restricting  eigenvalue  problems  to  those 
for  which  closed  form  solutions  exist. 

The  development  of  wave  propagation  theory  will  be  carried  to  the  point 
where  total  oblique  incidence  reflection  and  transmission  coefficients  can  be 
computed  as  a  function  of  frequency  for  rather  general  laminates  consisting  of 
laminae  of  the  treated  symmetries. 

The  mathematical  development  will  proceed  by  defining  the  complex  secular 
wave  equation  for  displacement  velocity  in  an  anisotropic  medium  and  trans¬ 
forming  it  to  allow  its  solution  in  terms  of  parameters  describing  wave  propa¬ 
gation  in  the  incident  medium.  The  solution  to  this  equation  provides  the 
complex  phase  velocity  and  displacement  velocity  eigenvectors  which,  when 
coupled  with  Snell's  Law,  yield  the  direction  of  the  propagating  waves. 

Using  this  information,  impedance  matrices  are  defined  which  relate  the 
stress  and  velocity  fields  for  the  propagating  waves. 

Once  the  stress  and  velocity  fields  are  characterized  for  the  media  on 
either  side  of  an  interface,  linear  boundary  value  equations  may  be  written 
relating  the  amplitudes  for  various  modes  of  wave  propagation  on  either  side  of 
the  interface. 

For  a  system  consisting  of  many  interfaces,  a  method  will  be  developed  for 
evaluating  the  total  reflection  and  transmission  coefficients.  This  technique 
is  a  generalization  of  one  developed  by  Scott1  for  normal  incidence  laminate 
propagation. 
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THEORY 

The  time  harmonic  equation  for  a  wave  in  an  anisotropic  medium  may  be  ex¬ 
pressed  in  the  form^. 

fc^iK  CKL  1Lj)vj  =  Py2  vi  •  t1) 

The  quantities  in  (1)  are  defined  as  follows: 

k  is  the  complex  wave  vector  for  the  propagating  wave 
u>  is  the  angular  frequency 

p  is  the  density 

is  the  vector  associated  with  the  i  th  component  of  wave  displacement 
velocity 

C^L  are  the  components  of  the  complex  elastic  constant  matrix  as  expressed 
L  in  contracted  engineering  notation  in  which  the  subscripts  run  over  6 
values,  and 

liR  is  a  matrix  of  the  form 


where  1;  are  complex  direction  cosine  defining  the  propagation  direction  of  the 
wave;  tne  transposed  matrix,  l^j,  is  defined  analogously. 

WAVE  EQUATION  IN 
AN  ORTHOTROPIC  LAMINATE 


For  the  most  general  case  to  be  considered,  an  anisotropic  solid  with  or¬ 
thorhombic  symmetry  and  principal  axes  along  the  x,  y,  and  z  directions  (1) 
takes  the  form: 
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where 


a  =  C^2.  cos2  0  cos2  <j>  +  cos2  0  sin2  $  +  C55  sin2  0 

8  =  C66  cos2  0  cos2  4  +  C22  cos2  9  sin2  4>  +  C44  sin2  6 

Y  =  C55  cos2  9  cos2  $  +  C44  cos2  0  sin2  $  +  C ^  sin2  9  (4) 

5  -  Cc12  ♦  C^)  cos2  0  cos  $  sin  <4 

e  =  (C13  +  Cgg}  cos  9  sin  9  cos  4 
£  =  (C23  +  C44)  cos  9  sin  0  sin  $ 

0  is  chosen  as  the  angle  of  incidence  of  the  wave  with  respect  to  the  plane 
of  the  laminate  (measured  up  from  the  plane)  and  4>  is  the  angle  between  the 
plane  of  incidence  and  the  first  principal  axis  in  the  plane  of  the  laminate. 

It  should  be  noted  that  a  number  of  possible  conventions  exist  for  choosing 
axes  in  an  anisotropic  medium.  In  a  unidirectional  fiber  reinforced  lamina  the 
first  coordinate  index  is  usually  chosen  along  the  fiber  axis  while  the  third 
index  is  chosen  perpendicular  to  the  plane  of  the  laminate.  On  the  other  hand, 
when  the  wave  equations  refer  to  a  crystallographic  system  the  third  axis  is 
conventionally  chosen  perpendicular  to  the  plane  of  hexagonal  symmetry  which 
would  be  the  counterpart  of  the  fiber  direction.  Hence,  when  values  of  elastic 
moduli  are  chosen  for  the  equations,  care  must  be  taken  to  ensure  that  they  are 
referenced  to  the  same  coordinate  system  as  the  equations.  In  this  paper  we 
adhere  to  former  convention. 

The  wave  equation  in  the  form  (3)  may  be  solved  in  closed  form  whenever 
the  angles  defining  the  propagation  direction  (9  and  <j>)  are  known.  This  is  not 
often  the  case  experimentally. 

In  most  practical  problems  involving  anisotropic  media  the  waves  present 
have  risen  from  refraction  or  reflection  at  a  given  interface  and  it  is  desired 
to  solve  the  equation  (3)  in  terms  of  the  parameters  of  the  incident  wave.  In 
order  to  treat  such  cases  the  complex  form  of  Snell's  Law 

k  cos  0  =  k'  cos  9'  (5) 

is  used,  where  in  general  9,  9',  k  and  k'  are  complex.  The  complex  angles  are 
necessary  foT  separately  defining  planes  of  constant  phase  and  attenuation 
which  may  not  always  be  parallel  in  problems  involving  oblique  incidence. 

By  substituting  the  relation  (5)  into  (4)  which  in  turn  is  substituted  in¬ 
to  (3)  a  secular  equation  is  derived  which  expresses  the  quantity  w2/h2 
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in  terms  of  the  known  quantities  6,  . ,  the  angle  of  incidence  9'  and  the  in¬ 
cident  wave  vector  k'.  ^ 


This  is  not  a  secular  equation  in  the  classical  sense,  since  the  off-diag¬ 
onal  elements  involve  the  variable  w/k;  however,  the  resulting  bicubic  equation 
does  provide  the  correct  eigenvalues  and  eigenvectors  for  the  refracted  waves. 

On  the  other  hand,  if  we  were  to  compute  actual  numerical  values  of  k  and  0 
from  (5)  their  substitution  into  (3)  and  (4)  would  yield  in  addition  to  one  cor¬ 
rect  solution,  two  spurious  solutions  corresponding  to  other  eigen  modes  trav¬ 
eling  parallel  to  the  true  refracted  wave. 


As  mentioned  above,  the  secular  equation  (3)  is  bicubic  in  n)2/k2  and  it  is 
easily  shown  that  it  is  also  bicubic  in  cos  0.  Expanding  (3)  we  have 


(°-£)  ((e-F1)  (”£) 

s[5Cv  )-5;] 

6  [55-*  (•-;?)] 


-  eC 


(6) 


From  the  relations  (4)  it  is  clear  that  cos  9  appears  everywhere  in  even 
powers  since  the  terms  in  cos  0  sin  0  (e  and  £)  always  appear  in  even  products 
and  in  fact  the  equation  can  readily  be  converted  to  a  cubic  equation  in  cos2  0 
using  the  substitution 


COS20 


cos 


2e' 


k'2 


(7) 


The  fact  that  the  equation  can  be  expressed  in  terms  of  even  powers  of 
cos  0  is  a  result  of  (3)  being  expressed  in  such  a  way  that  incident  and  re¬ 
flected  waves  move  with  the  same  speed. 


WAVE  EQUATION  IN  A  HEXAGONAL  MEDIUM 


For  the  case  of  hexagonal  symmetry  the  secular  equation  (1)  reduces  to  a 
biquadratic  equation  with  solutions  of  the  form 


(8) 
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1 


(10) 


These  are  denoted  as  the  pure  shear,  quasishear  and  quasi longitudinal  solu¬ 
tions  respectively,  where  a  is  the  angle  between  the  unique  principal  axis  and 
the  direction  of  wave  propagation.  This  solution  is  a  good  first  approximation 
for  the  problem  of  a  unidirectional  fiber  reinforced  laminate,  where  the  fiber 
direction  lies  along  the  principal  axis.  Clearly  this  solution  yields  isotropy 
in  the  plane  perpendicular  to  the  fiber  direction. 


For  the  case  in  which  a  wave  is  obliquely  incident  upon  a  hexagonal  medium 
(see  figure  1)  from  a  medium  of  known  velocity,  equations  (6),  (7)  and  (9)  can 
be  inverted  through  the  use  of  Snell's  Law  giving  for  the  pure  shear  case. 


^kY  _  (C66  '  C44^  cos2  e'  cos2  <t>  (k'2/w2)  +  p 

WPS  ^66  C11) 


_ Em _ 

k'2 

p  +  (C66  ‘  C44}  C°s2  *  C°s2  9' 

where 

k  is  the  magnitude  of  the  wave  vector  in  the  hexagonal  medium 

0  is  the  angle  of  refraction  in  the  hexagonal  medium 

p  is  the  density  of  the  hexagonal  medium 

9'  is  the  angle  of  incidence 
k'  is  the  wave  vector  in  the  incident  medium 

p  is  the  angle  between  the  x  axis  (fiber  direction)  and  the  plane  of 

incidence  (or  the  plane  of  refraction),  and  the 

Cji's  are  elements  of  the  elastic  constant  matrix  for  the  hexagonal  mate¬ 
rial. 
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For  the  cases  of  quasishear  and  quasi longitudinal  waves  the  expressions  are 
somewhat  more  cumbersome: 


cos2  9 


QS,  _  -B 

QL  7  + 


/!!  -  4  cn  c44\ 
\A2  A  ) 


(13) 


where 


and 


B  ■  (Cu  -  CJ4)  (Cu  -  2C44  .  C3J)  -  2(C13  .  C44)2 


COS^ 


A  = 


r 

r  2 

-  1 

C44  +  Clll 

“  P  „  +  (C11  ‘  C33)  C0S  * 

L  J 

[_k'  cos  9' 

J 

(cn  '  c33>  “s2  * 

2  cos2  9 

2 

_ 

- 

“1 

(C11  -2  C44 

+  C33)2  ‘  4(C13  +  C44> 

4 

COS  <j> 

’  E  D  A  N  C  E 

MATRICES 

AND  THE 

(14) 


(15) 


BOUNDARY  VALUE  PROBLEM 

In  order  to  solve  the  complete  boundary  value  problem  for  waves  traveling 
between  two  anisotropic  media  it  is  also  necessary  to  compute  the  normal  compo¬ 
nent  of  the  stress  fields,  or  traction  forces  associated  with  waves  on  either 
side  of  the  interface.  This  relation  is  given  by2 


where 


_T  =  niK  CKL  KLj 
in 

!i> 


(16) 


‘iK 


r  n  0  0 

x 


0  n 


0  n  n 
z  y 


n  0  n 
z  x 


0  n  n  n  0 
z  y  x 


(17) 


n.  are  normals  to  the  interface 
i 
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(18) 


Kj  are  wave  vectors  for  the  wave  within  the  given  medium  and  T.  (the  traction 
force)  consists  of  three  components  transforming  as  a  vector. 


The  quantity 


nlK  CKL  KLj 

u>  "  lj  (19) 

in  (16)  is  defined  as  the  normal  acoustic  impedance  matrix  and  is  useful  in 
solving  the  boundary  value  equations  for  a  given  interface. 


Boundary  value  equations  may  be  written  by  applying  the  conditions  that  the 
displacement  velocities  and  traction  forces  at  the  interface  between  two  media 
must  match.  This  implies: 


V  = 
u 


2  V 


u  u 


(20) 

(21) 


where  the  unprimed  variables  refer  to  the  incident  medium.  In  general,  for 
orthotropic  symmetry,  these  equations  will  contain  six  unknowns  corresponding 
to  three  possible  modes  for  the  reflected  and  refracted  waves  (2  quasishear, 

1  quasilon6itudinal) . 

Each  vector  may  be  written  in  the  form  =  Aj  vj  where  v^  is  an  eigen¬ 

vector  of  (3)  and  is  the  displacement  velocity  amplitude.  Since  the  compo¬ 
nents  of  each  v^  are  known  within  a  constant  factor  (10)  and  (11)  each  may  be 
broken  down  into  3  component  equations  allowing  solutions  for  the  six  unknowns 
to  be  obtained. 


Care  must  be  taken  in  (10)  and  (11)  to  insure  that  identical  components  are 
being  matched  for  the  incident  and  refracted  media.  For  media  in  which  princi¬ 
pal  axes  are  not  aligned,  this  is  most  easily  done  by  rotating  the  T's  and  v's 
on  the  right  side  of  the  equation  after  initially  computing  them  in  a  system  for 
which  the  coordinate  axes  and  principal  axes  are  coincident. 
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COMPUTATION  OF  THE  Z  MATRIX 

In  order  to  compute  the  Z  matrix  of  expression  (9)  it  is  necessary  to  know 
the  wave  vector  k  for  the  particular  mode  in  question.  For  refracted  waves  this 
information  comes  directly  from  the  solution  to  equation  (3).  When  reflected 
waves  are  considered  one  can,  by  the  use  of  symmetry  arguments,  obtain  the  k 
vectors  for  the  refracted  waves  directly  from  those  of  the  reflected  waves  by 
changing  the  sign  of  the  component  of  k  normal  to  the  interface.  If  all  three 
modes  of  the  incident  wave  are  present  this  will  permit  immediate  computation 
of  the  reflected  Z  matrices;  however,  if  only  one  mode  is  present  this  may  be 
used  to  determine  the  corresponding  reflected  mode  which  may  then  be  used  in 
place  of  the  incident  k*  and  9'  for  the  solution  of  equation  (3). 

Having  computed  the  oblique  incidence  reflection  and  transmission  coeffi¬ 
cients  for  waves  incident  upon  interfaces  between  anisotropic  media  of  interest, 
it  becomes  possible  to  compute  frequency  dependent  reflection  and  transmission 
coefficients  for  arbitrary  laminates  consisting  of  such  media. 

PLANE  WAVES  IN  LAYERED  MEDIA 


1 


2 


3 


n  +  1 


Figure  2.  Plane  Waves  being  Reflected  and  Mode  Converted  at  the 
Interfaces  of  a  Layered  Medium 


Figure  2  depicts  a  laminate  for  which  a  plane  wave  is  obliquely  incident 
upon  an  interface  from  the  lower  left.  The  wave  is  monochromatic,  and  by  defi¬ 
nition  infinite  in  width  and  may,  therefore,  be  equally  well  represented  by  any 
parallel  ray.  Upon  striking  the  interface  the  wave  produces,  in  the  most  gen¬ 
eral  case,  three  reflected  and  three  refracted  modes  whose  amplitudes  may  be 
denoted  by  the  vector 
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For  the  sake  of  convention  the  components  may  be  ordered  according  to  cor¬ 
responding  velocity  for  orthorhombic  symmetry.  For  hexagonal  symmetry  the  or¬ 
dering  quasi longitudinal,  quasishear,  pure  shear  is  less  ambiguous  since  modes 
do  not  change  at  points  where  velocities  cross  over.  It  is  clear  from  Snell's 
Law  and  elementary  geometry  that  for  plane  parallel  media,  even  numbered  multi¬ 
ple  reflections  in  a  given  direction  will  always  be  scattered  into  one  of  the 
three  modes  present  in  the  transmitted  wave.  In  other  words,  for  the  situation 
depicted  in  figure  2,  a  single  vector  a  as  in  (12)  is  adequate  for  describing 
all  of  the  waves  propagating  to  the  right  within  any  given  medium.  This  im¬ 
plies  that,  since  we  are  dealing  with  a  linear  system  of  plane  waves,  any  re¬ 
flection,  propagation,  transmission,  etc.  of  such  waves  can  be  treated  as  a 
matrix  multiplication  operation.  For  example,  a  vector  a  whose  components 
represent  the  three  possible  modes  of  propagation  with  a  given  (k *  sin  9')  for 
any  given  medium  will  have  components  of  the  form 

a^  =  ai  exp  i(wt  -  ki  •  r)  (23) 


Figure  3.  Waves  Referenced  to  the  Origin 

For  simplicity  the  phase  of  wave  (23)  will  be  referenced  to  point  P(r=0) 
(see  figure  3).  Thus  in  passing  from  medium  1  to  medium  2  the  vector  a  will 
undergo  the  component  transformation 

a*  j  =  Ttj  exp  i (wt  -  k'^  •  r) ,  (24) 

where 
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a  is  a  mode  vector  which  represents  waves  with  the  same  "k  sin  9"  inci¬ 
dent  at  P 

a^  is  the  complex  amplitude  of  the  i  th  component  of  such  waves  at  r=0 
and  t -0 

k^  is  the  wave  vector  for  the  i  th  component  and 

j  is  the  amplitude  transmission  coefficient  from  the  i  th  incident  mode 
to  the  i  th  refracted  mode 

The  variables  a',  k',  are  defined  analogously  for  medium  2  which  implies  the 
matrix  equation 

a'.  »  T-.  a-  or  a*  =  T  a  (25) 

j  1  ij  1 

It  should  be  noted  that  the  vanishing  of  the  phase  component  of  USlresults  from 
the  assumption  that  the  wave  front  travels  an  infinitesimal  distance  in  going 
from  medium  1  to  medium  2  in  the  neighborhood  of  P  and  is  independent  of  the 
fact  that  r  =  0.  Hence,  any  phase  shift  experienced  would  be  entirely  the  result 
of  transmission  phenomena  described  in  equations  (ICS)  and  (11)  from  which  Tjj  is 
computed. 

In  computing  phase  shift  matrices  corresponding  to  propagation  delays,  it 
is  convenient  to  describe  phase  shifts  seen  at  a  given  point  in  terms  of  fre¬ 
quency  domain  spectra.  Thus  a  wave  traveling  from  point  P  to  point  P'  in  medium 
2  will  undergo  a  phase  shift  such  that 

a'.  (P1)  =  a'.  (P)  exp  i(-k*  D  cos  9  '*) 

=  a^  (P)  exp  -  i(u>D  cos  9'^/Cj)  (26) 

which  is  described  by  a  diagonal  matrix  of  the  form 

$ij  =  5ij  exp  *i^wD  cos  9'i/ci)  (27) 

a'i  (f>,)  =  *ij  ai  (P)  (28) 

a  (P>)  =  7  a  (P)  (29) 

The  other  matrix  of  interest  is  that  associated  with  a  reflection.  Again 
no  propagation  delay  is  associated  with  this  phenomenon  and  reflection  matrices 
produce  the  transformation  in  medium  1 

(P)  “IXj  a.  (P)  (30) 

or 

b  =  R  a  where  b^  =  exp  i(oot  +  k^*r)  (31) 
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where  R;  is  the  amplitude  reflection  coefficient  for  i  mode  incident  waves 
going  toJ j  mode  reflected  waves. 

Again  the  coefficients  Rj;  are  computed  by  solving  the  equations  (10)  and 

(11). 

MULTIPLE  REFLECTIONS  IN  A  THICK  PLATE 


In  calculating  the  reflection  and  transmission  coefficients  of  a  thick 
plate  it  is  convenient  to  compute  the  total  amplitude  of  waves  at  the  points  P 
and  P'  and  to  consider  wave  fronts  as  being  reflected  between  those  points  (figure  3). 
This  is  permissible  since  it  is  easily  shown  that  the  phase  shift  of  the  wave- 
front  is  independent  of  the  point  from  which  the  reflection  is  assumed  to  occur. 

Then  for  an  incident  wave  A  (u)  the  total  reflected  wave  at  the  point  P  is  of 
the  form 


B(id) 

B(u) 

B(u) 


(R12  "  T21  *2  R23  *  T12  +  T21  *2  R23  *2  R21  *2  R23  *2  T12  ♦ 
+  *21  *2  R23  *2  R21  *2  R23  *2  R21  *2  R23  *2  *12  *  " A(uj) 


‘R12  *  T21  *2  R23  *2  So  «21  *2  *23  V"  T12>  A‘“> 

[*12  *  f21  *2  *23  *2  («-  521  *2  *23  ♦>''  fl2]  A(“> 


(32) 

(33) 


where  B(u>)  is  the  amplitude  vector  for  waves  of  frequency  u  traveling  in  medium 
1  with  a  negative  y  component  of  velocity  and  matrices  of  the  form  T. .  or  R.  . 
are  incident  from  medium  i  onto  medium  j .  }  } 


A  composite  reflection  coefficient  may  now  be  defined 


^R12  *  ?21  *2  R23  *2  '  R21  *2  R23  <t>2)  *12^ 


(34) 


which  describes  the  total  contribution  in  medium  1  from  interface  reflections 
occurring  to  the  right  of  medium  1  (see  figure  2) . 


In  a  similar  fashion,  the  wave  at  point  P'  emerging  into  medium  3  can  be 
written 


''M  -{T23  ‘2  T12 


+  T23  *2  R21  *2  R23 


‘12 


+  T23  *  R21  *2  R23  *2  R21  *2  R23  *  T12 


A(ui) 


(35) 
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*'<“>■  iT23  *  S  (R21  *2  *23  *2)l’  T12  }A<“> 

•  {  T2j  ♦  I  t  -  »21  *2  B23  ♦  I  ■'  f12  }*M  (36) 

The  total  transmission  coefficient  may  now  be  defined  as, 


REFLECTION  AND  TRANSMISSION 
COEFFICIENTS  FOR  LAYERED  MEDIA 


In  the  frequency  domain,  expressions  (34)  and  (37)  may  be  applied  itera¬ 
tively  to  construct  reflection  and  transmission  coefficients  for  multilayered 
media.  In  figure  2  a  wave  is  incident  from  the  left  upon  a  stratified  medium 
having  n  +  1  layers,  the  n  ♦  1  layer  being  the  incident  layer  and  the  n  th 
interface  being  the  first  interface. 


If  the  media  beyond  medium  n  and  in  front  of  medium  1  (with  respect  to  the 
direction  of  wave  motion)  are  considered  as_a  single  virtual  interface  having 
reflection  and  transmission  matrices  1L  Tn.j ;  then  the  problem  becomes  equiv 
alent  to  that  of  a  thick  plate  with  reflection  and  transmission  matrices, 


[r 


nf  +  Cnb 


t>  R  ,  (1 

n  n-l  n 


nb 


fr"  R  ,  4>  )_1  t  J 
n  n-l  nJ  nfJ 


(38) 


T 

n 


(I  - 


nf J 


(39) 


where  rnf  (rnb)  and  tnf  (tab)  are  the  forward  (backward)  reflection  and  trans 
mission  coefficients  respectively. 


Clearly  equations  (38)  and  (39)  form  recursion  relations  which  allow  the 
computation  of  and  Tn  in  terms  of  the  phase  shift  matrices  4>j  and  the  reflec 
tion  and  transmission  matrices  t\f,  tjb,  r^f,  r^b. 


This  set  of  recursion  relations  can  be  put  in  matrix  form  using  the  method 
of  Scott1  by  employing  the  definitions 


„  -l 

R_  »  a  b  , 
n  n  n 


a  */rrt~-r-trr.$  R  ,  $  ♦  t  .#  R  ,  <j>  ^ 

n  I  nf  nf  nf  nf  nb  n  n-l  n  nb  n  n-l  nl 


(40) 

(41) 


and 


bn  *  ^nf  ‘  cnf  rnb  *n  ^n-i  *n^ 


(42) 
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where  the  bar  matrix  notation  is  suppressed.  Expression  (40)  may  be  shown  to 
follow  from  (38),  (41)  and  (42)  by  direct  substitution.  By  use  of  definition 
(40)  we  may  write 


R 

n 


a 

n 


x 


*nf  ' 


r  f  t  i  r  ,  $  R  <{>  )$  1  b  +  t  . 
nf  nf  nb  n  n-i  n  n  n-l  nb 


-1 


=  a  1  b  ,  b  1 ,  <J>  b~ 1 
n  n  n-l  n-l  n  n 


(43) 


Thus  for  each  pair  of  matrices  a^  bn  a  matrix  relationship  may  be  established 
of  the  form: 


^nb  *  rnf  cnf  rnb>  V  rnf  cnf  ♦n1  \ /Cl 


_tnf  rnb  ^n 


nf  vn 


(44) 


Extending  the  relationship  through  all  of  the  layers  of  the  laminate  yields: 

.  l 

r  .  1  4  r_r.d>-l 


^mb  "  rmf  tmf  rmb^  %  rmf  tmf  ^in1 


m=n ,  T 


"''rnf  rmb  V 


-1  -1 
tmf  ^m 


(45) 


The  expressions  (44)  and  (45)  above  could  be  expanded  into  a  form  contain¬ 
ing  6x6  and  6x3  matrices;  however,  the  above  forms  are  much  easier  to  work  with 
and  are  expressed  in  terms  of  the  3x3  matrices  desired  in  the  result. 

From  expression  (39)  a  similar  definition  can  be  made  for  transmission 
matrices.  Thus,  letting 


c 

n 


"n-l  yn 


(46) 


and  as  before 


b  * 
n 


t'i  r  ,  ip  R  ,  $  ) 
nf  nb  n  n-l  n 


the  expression  (39)  becomes 


T 

n 


This  leads  to  the  recursion  relation 


(47) 


(48) 


13 


NADC- 80099-60 


T  =  [c  b'\  $  1  [t‘*  -  t*i  r  .  *  a  .  b'1.  *  l*1 

n  1  n-1  n-1  nJ  1  nf  nf  nb  n  n-1  n-1  nJ 


=  cn_l  '  [t*'  b 


nf  rn  “n-1  *  *nf  rnb  ^n  an-l^ 


-1 


from  which  is  derived  the  matrix  relation 


(t„u  -  r  .)$_  r_x.  t 


0  -t'i  r 


nb  *nf  nf  nbJ*n  nf  nf  ^n 
nf  *nb  *n 


(49) 


(50) 


A  matrix  product  of  the  type  in  expression  (45) could  readily  be  constructed,  how¬ 
ever  it  is  clear  that  there  is  no  mixing  of  the  c  s  with  other  matrix  elements 
and  the  immediate  result  is 


T™  •  <if  •  c 


(51) 


where  b^  is  obtained  directly  from  expression  (45) . 


RESULTS 

The  above  analytical  technique  has  been  utilized  for  calculating  reflection 
and  transmission  coefficients  for  a  number  of  layered  anisotropic  and  isotropic 
materials  both  lossless  and  attenuating.  Although  the  author  has  not  as  yet 
found  independent  results  in  the  literature  for  checking  the  most  general  cases 
which  can  be  treated  by  this  technique,  verification  has  been  performed  for  a 
number  of  results  for  reflection  and  transmission  between  crystalline  media  and 
crystalline  and  amorphous  media  given  by  Auld.2  in  addition,  formulas  given  by 
Ewing  Jardetzky  and  Press-*  for  liquid  to  solid  and  solid  to  liquid  interfaces 
have  also  been  checked  and  give  agreement  of  the  type  shown  in  figures  4  and  5. 
Excellent  agreement  was  also  found  with  results  given  by  Brekhowskikh4  for  trans 
mission  through  a  thick  aluminum  plate  in  water  as  a  function  of  angle.  Further 
more  all  solutions  generated  for  very  general  laminate  systems  have  satisfied 
energy  conservation  laws. 


CONCLUSIONS 


The  above  formalism  is  believed  to  provide  an  accurate  and  efficient  tech¬ 
nique  for  treating  rather  general  problems  involving  the  oblique  indicence  trans 
mission  and  reflection  of  stress  waves  through  anisotropic,  layered  media.  The 
technique  has  proved  adequate  to  handle  problems  involving  evanescent  waves  and 
material  attenuation  and  should  be  useful  in  treating  practical  problems  involv¬ 
ing  composite  materials. 


REFLECTION  COEFFICIENTS 
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Figure  5  -  Reflection  Coefficients  for  Waves  Traveling  from 
a  Solid  into  a  Liquid 
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